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Abstract 

The Einstein-scalar field theory can be used to model gravitational 
physics with scalar field matter sources. We discuss the initial value 
formulation of this field theory, and show that the ideas of Leray 
can be used to show that the Einstein-scalar field system of partial 
differential equations is well-posed as an evolutionary system. We 
also show that one can generate solutions of the Einstein-scalar field 
constraint equations using conformal methods. 



1 Introduction 

In Newtonian theory, one models gravitational physics by studying a lin- 
ear elliptic Poisson equation for the Newtonian potential on a fixed absolute 
background space and time, with the motion of material bodies governed 
by the Newtonian force equation on this fixed background. By contrast, in 
general relativity the gravitational field is modeled using Lorentzian space- 
times whose curvature reflects the material and field content of the spacetime. 

*Partially supported by the NSF under Grant PHY-0354659 
tPartially supported by the NSF under Grant DMS-0305048 



1 



Mathematically, a Lorentzian spacetime is a pair (M n+1 , g) where M n+l is a 
smooth manifold of dimension n+1 (in everyday physics n — 3, but higher di- 
mensions are sometimes considered for modelling electromagnetism and other 
interactions with gravity), and g is a pseudo-Riemannian metric of signature 
( — , +,..., +). The metric distinguishes timelike directions g(X,X) < for 
a tangent vector X (along the possible path for a massive physical object), 
null directions g(X,X) = (along the possible path for a massless physical 
particle), and spacelike directions g(X,X) < 0. The physical time between 
a pair of events, as would be marked by a proper clock, corresponds to the 
g-length of the timelike trajectory followed by that clock. The path followed 
by a test particle which is free of non gravitational forces corresponds to a 
timelike geodesic in the spacetime. 

It is well known that theoretical studies of general relativity have pre- 
dicted such strange and interesting phenomena as the expansion of the uni- 
verse, black holes, gravitational lenses, and gravitational waves. These phe- 
nomena, all of which have now been confirmed either by direct or indirect 
observation, were originally discovered via studies of solutions of the Ein- 
stein field equations. It is also well known, at least among mathematicians, 
that Einstein's equations present a number of very challenging mathematical 
problems, such as the cosmic censorship conjectures and the question of the 
nonlinear stability of black holes. 

In this brief review, we show that in two important ways, the ideas of 
Leray have played an important role in the study of Einstein's equations. In 
previous studies these ideas were applied in the absence of a scalar field. Here, 
since scalar fields are now viewed as possibly important for understanding the 
apparent acceleration of the expansion of the universe, and since including 
them results in some additional interesting features in the analysis, we work 
with the Einstein-scalar field system. We introduce this system in section 2 
and discuss its Cauchy formulation in section 3. We then show in section 4 
how Leray's ideas play a role in understanding the evolutionary aspect of the 
Einstein-scalar field system. In section 5 we prove the existence of solutions 
to the Einstein-scalar constraint equations via the conformal methods, where 
seminal ideas of Leray regarding solutions of nonlinear elliptic equations have 
played a crucial role. 

We end this section with a note about notation. Wherever there may arise 
a possible confusion, we use bold faced symbols for spacetime variables and 
tensors. For initial data sets, we reserve the notation of over-barred symbols 
for physical variables, which satisfy the relevant constraint equations, and 
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may or may not be "time-dependent" depending on the context. Unadorned 
symbols are used for free conformal data as described in section 5. Finally 
the occasional appearance of a "tilded" -symbol (as in (15 .4p below) appears 
when we need to introduce an intermediate quantity which is neither a free 
variable, nor a physical one. 

Acknowledgement: The authors would like to thank the Isaac Newton 
Institute of Mathematical Sciences in Cambridge, England for providing an 
excellent research environment during the program on Global Problems in 
Mathematical Relativity during Autumn, 2005, and again in October 2006, 
where some of this research was carried out. 



2 Einstein - scalar field equations. 

For general source fields, the Einstein gravitational field equations take the 
tensorial form 

G(g) = T($,g), (2.1) 

where G(g) is the Einstein tensor, which is a second order differential oper- 
ator on the metric defined by G(g) := Ric(g) — |Rg with Ric(g) denoting 
the Ricci tensor of g and R denoting the scalar curvature of g, and where 
T($, g) is the stress-energy or energy-momentum tensoi0, a specified func- 
tional of the matter source fields $ and the metric. The specific form that the 
stress-energy tensor takes depends upon the matter source fields presumed 
to be present in the physical system being modeled. Here, we presume that 
$ is a scalar field, which we label with potential function V(ty), and we 
set 



(2.2) 



In addition to the equation (12. ip . the Einstein-scalar field theory includes 
a field equation for which reads 

dV 

V a d a V = — , (2.3) 

where V denotes the covariant derivative compatible with g. While this extra 
equation may simply be added to the theory by hand, it also follows directly 



1 Note that we have chosen units so that 8ir times the Newtonian gravitational constant 
is set equal to one. 
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as a necessary consequence of equation (12.11) together with the geometric 
properties of the Einstein tensor: One readily verifies that the Bianchi iden- 
tities for the curvature imply that the Einstein tensor satisfies the identity 

div g G(g) = (2.4) 

where div g denotes the divergence operator for the metric g. This condition 
together with (12. ip implies the conservation law 

div g T = 0. (2.5) 

We readily verify that (I2.5P applied to the scalar field stress-energy tensor 
(12.21) results in the field equation (12. 3p . 

How does one choose the scalar field potential function V(^)7 While 
there are many possibilities, we note that V(^) = y\I/ 2 corresponds to the 
massive Klein-Gordon field, while setting V(^) = A for a non zero constant 
A and requiring that \1/ = produces the vacuum Einstein theory with non 
zero cosmological constant A. 

3 The Cauchy problem: constraints and evo- 
lution. 

The Einstein-scalar field system of partial differential equations on a 3 + 1 
spacetime consists of eleven equations (12. ip and (12.31) for the eleven field vari- 
ables g^ and \1/ (for n+ 1 dimensions, there are h(n+ l)(n + 2) + 1 equations 
for the same number of field variables). One of its most characteristic fea- 
tures, however, is that it is both an under-determined and an over-determined 
system, in the sense that if one formulates the Einstein-scalar-field system 
as a Cauchy problem, there are constraint equations which must be satisfied 
by any candidate set of initial data, and as the data evolves there are certain 
of the field variables whose evolution is entirely at one's discretion. Both of 
these features reflect the spacetime covariance of the theory (i.e., the theory 
has the spacetime diffeomorphism group as its gauge group). 

To see these features explicitly, we now sketch out ann + 1 decomposi- 
tion of the Einstein-scalar field variables and equations. Given a spacetime 
(M n+1 , g), we start by choosing ann + 1 foliation of the spacetime manifold 
F t : S n -> M n+1 (t G M), for which each of the leaves F t (E n ) of the foliation 
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is presumed spacelike. We also choose a threading of the spacetime by a con- 
gruence of timelike observer paths T p : R — > M n+1 (p e S n ). The choice of 
a foliation and a threading, together with a choice of coordinate patches for 
£ n , automatically determines local coordinates (x° — t, x 1 , . . . , x n ) and local 
coordinate bases (J^, . . . , covering M n+1 . We may then, without loss 
of generality, write the metric (locally) in the form 

g = -N 2 9 l ®9 l + <g> 9 j 

where (9 t = dt, 9^ = dx^ + ftdt) is the one-form basis dual to the surface- 
compatible tangent vector basis (e± = -^(^ - {3 j -£j), gfr • • • , gf^) with e± 
the normal vector field to F t (T, n ). Here N is the positive definite "lapse 
function" , ft are the components of the spacelike "shift vector" , and 7^ are 
the components of the spatial metric tensor. We note that for each choice 
of t, 7(i) = jij(t)dx l <8> dxi is the induced Riemannian metric on the leaf 
F t (S n ). We use the notation K{t) = Kij(t)dx l <g> dx^ to denote the second 
fundamental form defined by the foliation. 

For the scalar field there is no need to do any space + time decom- 
position. However, we shall use the notation rp to denote the restriction of 
\I> to one of the leaves of the chosen foliation, and we shall use the definition 
^ := iHJ^ ~ ^dx™^ ^ or convemence m working with the time derivative 

If we now use the usual n+ 1 decomposition to express the spacetime cur- 
vature in terms of the time dependent spatially covariant quantities 7, K, N, (3 
and their various derivatives and (spatial) curvature, we find that the Einstein- 
scalar field equations (12.1 j) and (12.31) split into two types: constraint equations 
which require any choice of initial data to satisfy certain identities, and evo- 
lution equations which describe how the spatial fields evolve from one leaf of 
the foliation to the others. Explicitly, we have the following: 

Constraint Equations 

From the G±± equation derived from (12.11) . we obtain the Hamiltonian con- 
straint 

2N~ 2 G L± = Ry- \K\* + {trKf = tt 2 + |V^| 2 + 2V($). (3.1) 

From the G±j equations derived from (I2.ip . we obtain the momentum con- 
straint 

- N- 1 G ±j = V m Kf - djtrK = TfOjtp. (3.2) 
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Note that these equations constrain the choice of the data (7, K, ip, n); they 
do not involve the lapse and shift. We refer the interested reader to [2] for a 
survey on the constraint equations. 

Evolution Equations 

From the Gy equations derived from (12. 11) . we obtain 

^Kij = N(R ij -2K im K™ + trKK ij -dSd j $ + : ^%V$)) 

-LvfiN + CpKij, (3.3) 

where here and above V is the covariant derivative associated to 7, Rij are 
the components of the spatial Ricci tensor calculated from 7, and C denotes 
the Lie derivative operator. This is an evolution equation for K. We obtain 
an evolution equation for 7f from the spacetime field equation (12.31) for t/j 

r\ it t 

-7f = N(A^ + tr Ktt - -j) + j mn d m Nd n i; + CpTt, (3.4) 
01 dip 

where A^ is the Laplace-Beltrami operator for the metric 7. 

We have evolution equations for K and 7f. What about 7, N, (3 and ipl 
The evolution equation for i/> comes from the definition for 7f: 

^ = Nn + Cfrf>. (3.5) 

The evolution equation for 7 comes from the definition of the second funda- 
mental form: 

^ tJ = ^NKij + C^j. (3.6) 

For the other field variables, iV and j3, there are no evolution equations. 
This freedom to choose iV and (3 and their evolution any way one wishes 
reflects the gauge invariance of the field equations under the action of the 
diffeomorphism group. 

To summarize, the Cauchy formulation of the Einstein-scalar field equa- 
tions asks that one choose the initial data set (7, K, ^},tt) subject to the 
constraint equations (I3.ip - fl3.2l) . One then chooses N and (3 freely in time, 
to fix the gauge, and finally one proceeds to evolve (7, K, ijj, 7f) via the evo- 
lution equations just listed. Note that if one chooses the constraints to hold 
initially, in any accurate evolution they must remain satisfied for all time. 
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4 Evolution system. 



Does an initial data set satisfying the constraint equations (I3.1l) - (l3.2p always 
generate (via the evolution equations (I3.3p - (I3.6I) ) a spacetime solution of 
the Einstein-scalar field equations (12. lft - ( 12731 ? To show this, one needs to 
prove that the system is well-posed in some appropriate sense. Here we 
discuss the well-posedness results of Leray-Ohya and of Leray, as applied to 
the Einstein-scalar field system. Note that our focus here is solely on the 
evolution system. As in the vacuum case for the wave gauge [10], one must 
use the Bianchi identity to show that the constraints are preserved in the 
evolution, thus yielding a local existence result for the full field equations. 
We refer the reader to [1] where this issue (again for vacuum solutions) is 
addressed in a similar setting to the one considered here. 

To check if a given system is hyperbolic in the Leray-Ohya sense, one 
seeks to diagonalize the matrix of the principal parts (highest derivatives) 
of the system. If this diagonalization can be done, the system is called 
"quasi-diagonal" , and is a Leray-Ohya hyperbolic system if in addition each 
operator has a characteristic cone which contains the metric cone. As a con- 
sequence, it can be shown that it is well-posed in Gevrey classes of functions. 
These are spaces of C°° functions whose successive derivatives satisfy certain 
inequalities which are generally too weak to imply the convergence of the 
corresponding Taylor series. Well-posedness means that a set of initial data 
in such classes will generate a spacetime solution with evolving data in the 
same classes. It also implies causal propagation, with the domain of depen- 
dence of the data determined by the causal cones of the spacetime metric, 
as well as continuous dependence of the evolved solution on the choice of 
data. Note that the principal parts of a Leray-Ohya system may generally 
have multiple characteristics. Also note that in verifying the criteria for a 
Leray-Ohya system, one need not have the same order for each of the various 
evolution equations which make up the system. 

If in fact the operators in the principal part of a Leray-Ohya system do 
not have multiple characteristics - i.e., if there exists a cone in the cotangent 
plane such that each straight line passing through a point in its interior cuts 
the characteristic cone in iV distinct points if the operator is of order N, then 
the system is called Leray hyperbolic. It can be shown that such a system is 
well-posed in Sobolev spaces as well in Gevrey class spaces. 

For a more detailed discussion of Leray-Ohya and Leray hyperbolicity, 
see H31 HH S]- 
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We now apply these ideas to the Einstein-scalar field system. In doing 
so, we work with the Einstein-scalar field system in mixed first order (for 7 
and K) and second order (for x/j), to mesh with the extant treatments. 

4.1 Leray - Ohya hyperbolic system for */,K and ip. 
We set Bq := — Cp and we consider the system 

Bo Rij — B ViRjo — B VjR i0 + VjVj-Roo — Fy 

which is derived by taking linear combinations of equations from the Einstein- 
scalar field system (see pQ for a similar derivation in the vacuum case). If 
B jij is replaced by its value in terms of K, 

d % = -2NK tv (4.1) 

this system reads as a third order system for K of the form 

B (-N- 2 B 2 Q + V h V h )K i:j = fa (2 in 7, 2 in K, 3 in N, 3 in 4>) + f tJ , (4.2) 

with 

f ir .= N\7 t d 3 (N- 2 Bl-W h d h )N. 

The numbers appearing above with (and below with h) denote the order 
of the highest derivatives of the unknowns which occur in that term. The 
additional term NVidj(N~ 2 BQ — V dh)N is clearly third order in 7 (and 
fourth order in N). 

The wave equation for ip reads, in terms of 7, K and the presumably 
specified variables iV and (3, as follows: 

- - - dV 

-N- l d {N- l d ij) + N- l f j Vi(Nd^) + KtN^d^ = - T . 

dip 

Applying the operator d and using (14.11) gives an equation of the form 

d (-N- 2 d 2 + V h d h )i) = h (1 in 7, 1 in K, 2 in N, 2 in (4.3) 

The principal matrix of a system of partial differential equations Eb{ua) = 
0, with unknowns ua, is obtained by assigning to each unknown an integer 
m(uA) and to each equation an integer n(E B ) such that the highest deriva- 
tives of ua appearing in Eb are at most of order m(it^) — ti(Eb)- The 
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principal part relative to ua in the equation Eb = then consists of the 
terms of order m{uA) — n{Es) in ua- It is zero if there are only terms of 
smaller order. These integers m(uA) and ti(Eb) are collectively called the 
Leray-Volevic indices for the system. 

For the system of PDEs (TO) and with the field variables 7(f), 

and ?/>(£), we choose the Leray-Volevic indices 

m (fy) = 3, m(K) = 3, m(#) = 4, (4.4) 

n(TCID = 2, = 0. nP3J = l. (4.5) 

The principal matrix is then a triangular matrix, with the elements in the 
diagonal being the derivative Bq along the normal e± to the space sections 
and the product of this operator by the wave operator in the spacetime 
matrix. These operators are causal and hyperbolic. Both are such that their 
characteristic cones at a point contain the metric null cone. However the 
non-diagonal form of the principal matrix does not permit one to conclude 
that it is Leray hyperbolic. 

If we replace K in (14.21) by its value in terms of 7 using (14.11) . and give to 
7 the index 4, we obtain for (14.21) and (14.31) a diagonal system with principal 
operators the wave operator and the operator 

dl(-N- 2 8* + V h V h ). 

This operator has a double characteristic, the spacelike hyperplane, and the 
system is only Leray-Ohya hyperbolic, relative to the Gevrey class of index 2. 
The Cauchy problem for this system is well posed in this class; the domain of 
dependence of the solution is determined by the light cone of the spacetime 
metric. 

4.2 Leray hyperbolic system for 7, K and ip with a 
lapse condition. 

The system for K and ip can be put into Leray hyperbolic form if we 
impose a condition on the lapse function N which makes it a quasi-diagonal 
system for 7, K, $ and now also N. 

To remove the term /y, which introduces non diagonal elements into the 
principal matrix, we require that iV satisfy a wave equation with source term, 
with that source term being an arbitrarily specified function F 

-N- 2 d* N + V%N = F, 
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and we then insert this equation into (I4.2p . We now consider the system 
consisting of (14. 1 1) - 04. 3 p with this change, together with the equation 



obtained by taking the do derivative of (14.21) . and using (14. lft . We choose for 
N and (14 .6p the Leray-Volevic indices 



The system is now quasi-diagonal with hyperbolic diagonal elements given 
by do(—N~ 2 dl + V 4 V;) and B . This leads to the following result. 

Theorem 4.1 The system fl^. i| ) ; {^.Ity , and flV-fip is a Leray causal 

hyperbolic system for 7, K, if) and N. 

5 The constraints 

The Einstein-scalar field constraints consist of the n+1 equations (I3.1l) - (l3.2p . 
to be satisfied by the initial data (7, K, x/j, 7f) on an n-dimensional manifold 
E. Locally, due to the symmetry of these tensors, this initial data can be 
regarded as a set of n(n+ 1) +2 functions, which makes the underdetermined 
nature of the constraint equations apparent. We recall that (7, K, V>,7f) de- 
notes a set of initial data which satisfies the constraint equations; we use the 
same quantities without the over bars to denote functions which we choose 
freely in order to construct the data (7, K, ■0, 7f). Also for convenience here, 
we restrict our considerations to the so called cosmological case, for which E 
is a compact manifold (see [7] for a treatment of the asymptotically flat case). 
Even in vacuum there are infinitely many solutions of the constraints, de- 
pending on arbitrary transverse-traceless (divergence and trace free) tensors, 
which can be interpreted as "radiation data". 

5.1 The conformally formulated constraints 

The conformal method involves decomposing the data (7, K, tt) into cer- 
tain parts which are chosen freely, and other parts which are determined by 
solving equations which we derive from the constraint equations. We consider 
the case n > 3. The most basic piece of the freely chosen data is a choice of 



d ((-N- 2 d 2 +V l d i )N) = d F, 



(4.6) 



m(iV)=4, rifTO = 1. 



(4.7) 



10 



a Riemannian metric 7, or rather the conformal class of metrics represented 
by 7. The physical metric 7 is required to be conformally related to 7. One 

sets 

4 

7 = (Yjn-2 7 

for a positive function 1/7 on S. The following identity then holds between 
the scalar curvatures of 7 and 7 



n + 2 



4(71 — 1 ) 

i?(7) = -<p-iZ$ ( _A__2a 7 ^ - i?( 7 V ) . (5.1 ) 



On the other hand, the divergences of traceless contravariant symmetric 2- 
tensors are related by the identity 

2(n + 2) 

div 7 P = »-a div 7 .P (5.2) 

if 

_ 2(n + 2) 

P = (p n-a P. 

Using (15.21) (applied to the traceless part of K) together with (15.11) . the 
Einstein-scalar field constraints may be conformally reformulated as follows. 



5.1.1 The Momentum constraint 

The momentum constraint may be expressed with respect to the background 
metric ip by 

n — 1 _2n_ 2(n + 2) _ 

div 7 fT = <£>™-2Vt + <£> ™- 2 J, (5.3) 

n 

where r = ti = tr^K is the mean curvature, the physical extrinsic curva- 
ture (second fundamental form) K is related to K (as contravariant tensors) 
by 

2(n + 2) ~ J 1 

K — ip ™- 2 K + -7 
n 

with 7 _1 denoting the contravariant form of the metric 7, and where J : = 
— ttV^. We have shown in [7J, following ideas originating from York [18], that 
it is useful to associate to the background conformal metric 7, a function iV 
which is related to the original lapse by the equation^] 

2 This relation consists in requiring that each metric has as associated to it an initial 
lapse with the same "densitized lapse" . 
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or simply 

2re ~ 
N = if~N. 

Then the "physical" scalar field initial data (t/>,7f) consists of ■?/> = ip and 

7f = rN~ 1 doip = (p~~^7i where n = N^doip, 
and it follows that 

_ 2(n + 2) 2(n + 2) 

J = — (f "- 2 7rV"0 = ^ "- 2 J, where J = — wVip. 

Hence the Einstein-scalar field momentum constraint equation does not con- 
tain (p if Vr = 0. Setting 

(5.4) 

with C l con fX the conformal Lie derivative of 7 (or conformal Killing opera- 
tor) with respect to a vector field X and with U a freely specified traceless 
2-tensor, the system may be regarded as a self-adjoint linear elliptic system 
for X, as follows. 

n — 1 22_ 

A 7 conf X = -div^/7 H ip~ Vr - ttV^. (5.5) 

n 

In summary, if we begin with a choice of "free" initial data (7, U, r, ip, n) and 
solve the conformally formulated momentum constraint equation (15.51) to 
determine K as indicated in fl5.4[) then K satisfies the momentum constraint 
equation (15. 3p . 

5.1.2 The Hamiltonian constraint 

If we specify the initial data (7, U, r, ip, n) and use the identities (I5.2p and 
(15.11) . then the Hamiltonian constraint equation ( 13. ip becomes a semilinear 
elliptic equation, called the Lichnerowicz equation, for ip. We have shown 
[H [8] that it takes the form 

H = A^tp-f(<p) = 0, (5.6) 

with 

. . 3n-2 n + 2 
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where we set a 



n-2 



and we let 



4(n-l) ' 



IWI?) , 




) 



and 




■ii 



( 



n — 1 



n 



t 2 - 4y(^) 



We observe that ^ w > 0, while the sign of B T ^ depends on the relative 
values of r and V(i/j). Note that in the constant mean curvature case the 
system of equations are "semi-decoupled" in that we may first solve the 
momentum constraint equation (15.31) (which does not involve ip) and then use 
the resulting K to formulate the Lichnerowicz equation as described above. 
If we can find a positive solution ip to this equation then this determines the 
physical metric 7 and second fundamental form K as well as the physical 
initial data 7f) for the scalar field. 

In [3 E] we establish a number of results regarding the existence or non 
existence of solutions for the system consisting of (15.51) and (15.61) . We present 
here some of the existence results, in a low regularity setting, on manifolds 
with r constant (constant mean curvature, or "CMC", initial data). The as- 
sertion that these results hold in a low regularity setting follows from meth- 
ods established by Choquet-Bruhat [3] and Maxwell [151 EE] for the vacuum 
Einstein constraint equations. 

5.2 Existence theorems. 

We denote by Wf the usual Sobolev space on (£, 7), consisting of functions 
with (for s a positive integer) all weak derivative of orders less than or equal 
to s lying in L p , and by Mf the space of Wf Riemannian metrics (which 
consists of an open cone in the space of all Wf 2-tensors if s > 2 and p > |, 
or if p = 2 and s > |). We denote W 2 S by H s and M 2 S by M s . 

5.2.1 Solving the momentum constraint. 

Given a traceless tensor U on E, the conformally formulated momentum 
constraint equation (15.51) is a linear elliptic equation for the vector field X. 
The kernel of A 7jCon j consists of the space of conformal Killing vector fields. 
The inhomogeneous term div 7 £7 is orthogonal to this space. The proof of the 
following theorem follows from known theorems for linear elliptic systems. 
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We suppose that s is an integer We let F denote the left hand side of (15.51) 
in the CMC setting, so that 

F = -div 7 U - TrVvf). 

Theorem 5.1 The conformally formulated momentum constraint equation 
< f5T5)) with 7 G M\, p> § has a solution X G W$, 1 < q < p, if F G L q and 
if J = 7T is orthogonal in the L? sense to the space of conformal Killing 
(CK) vector fields on 

Moreover the solution is uniquely determined up to the addition of a con- 
formal Killing vector field. There exists a constant c(j), depending only on 7, 
such that the unique solution which is orthogonal to the space of CK vectors 
satisfies 

\\X\\ w ,<c^)\\F\\ Lq . (5.7) 
Corollary 5.2 // 7 G M s 2 , s > § and F G H s _ 2 , then X G H s . 

Proof. We first remark that if s > ~, then the Sobolev embedding 
theorem C H s for p < w _^" +4 implies that if s > ~ there exists p > | 
such that the embedding holds. 

The corollary is then established in the usual way, by differentiating the 
equation and using the Sobolev multiplication and interpolation properties. 



5.2.2 Solving the Hamiltonian constraint. 

Satisfying the Hamiltonian constraint is equivalent to finding a positive so- 
lution of the Lichnerowicz equation (15.61) . To prove the existence of positive 
solutions we use the method of sub and supersolutions. The early approaches 
to solving the Lichnerowicz equation made use of the Leray-Schauder de- 
gree. The method employed here uses estimates for linear elliptic equations 
together with the Arzela-Ascoli theorem; see for example [TT]. The first so- 
lutions were found in Holder spaces, then in Sobolev spaces H s , s > + 1. 

3 The case of non integer s > § is treated by Maxwell [16] for the vacuum Einstein 
constraint equations. This of course requires working with distributional solutions if n=3. 
For integral choices of s > Q we in particular have s > 2 when n > 3. Thus the formulas 
presented here involve pointwise almost everywhere defined derivatives. 
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The regularity has since been reduced to W$ , p > § [3 J and to H s , s > § 
[T6] . in the absence of scalar field. 

We give a general theorem, adapted to equations of the type of the Lich- 
nerowicz equation, with or without a scalar field. Consider the semi-linear 
equation 

A 1 <p = f(x,<p)= Q i( x )^' ( 5 - 8 ) 

i=l,...,N 

on the compact manifold (E, 7), where x G E, and pi G R. We say that 93 _ is 
a subsolution of (15. 8p if A 7 </?_ > /(x, </?_) and y? + is a supersolution of (15.81) 
if A 7 (^ + < f(x,(p+). 

Theorem 5.3 Equation / 1 5. 8\) admits a positive solution tp G W%, P > § , 
provided the following conditions are satisfied 

(a) 7 G Mf ; p > § , and a { e V for i = 1, . . . , N. 

(b) The equation admits a strictly positive subsolution and supersolution 
ip + , both in W%, with < < y? + < 00. 

T/ie solution if then satisfies y?_ < < and zs unique if f(x,y) is 
monotonically increasing in y for each x G E. On t/ie ot/ier hand, if all of 
the ai are of the same sign then no positive solution exists. 

Corollary 5.4 // in addition 7 G M s , s > |, 04 G i/ s _2, i=l,...,N, then 
V? G i/ s . 

It can be proved by using conformal invariance that, in the case of the Lich- 
nerowicz equation, the uniqueness of the solution is independent of the sign 
of (see Theorem 7.12 of [1]). 

In the original analysis of the Lichnerowicz equation for vacuum CMC 
data on compact manifolds, the full Yamabe theorem^ is employed to fix the 
sign of the linear zero order term [11] . A verification that one only needs to 
control the sign of the scalar curvature (a much easier result), even in the 
low regularity setting, is provided by results of [5] and [IB] . The following 
result provides the analog of this control in the presence of a scalar field (see 
Proposition 1 of [8]). 

4 This says that every metric on a compact manifold is conformal to one with constant 
scalar curvature. The proof of this theorem was completed by Schoen [17] after essential 
contributions by Yamabe, Trudinger and Aubin. We refer the interested reader to [TJl E] 
and the references contained therein. 
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Theorem 5.5 (The Yamabe-scalar field conformal invariant) The func- 
tional on Hx(E) (for 7 G Mf and tj) G given) defined by 




admits an infimum, 3^ ([7]) > —00, which is a conformal invariant. Its sign 
determines the Yamabe-scalar field classes of pairs (7, VO- A pair { r y ) ip) with 
^([7]) < (respectively 3^ ([7]) = 0, or 3^ ([7]) > 0) can be conformally 
transformed to a pair such thatH^ < (respectively IZy^ = 0, orlZ^^ > 0) 
on E ; and moreover if IZ^^ maintains a fixed sign on E it is necessarily of 
the same sign at ^([7]). 

The proof of the following existence theorem relies on the construction of 
sub and supersolutions tp_ and ip + . We have supposed that r is a constant 
and that V is a smooth function of i/j G Wf C C°(E), since p > |. We 
therefore also have B T ^ G C°(E) C 

The expression ^4 7 ^ = c n (\K\^ + 11 2 j shows that A y % n > 0, and that, 
for a solution K of the momentum constraint, -A x n G L p since 

Wf x Wf c L p when p > -. 

2 

We present results here in the case that B T ^ = c n (^^-t 2 — 4V( , 0)) > 0. We 
refer the interested reader to [S] for a more general treatment. 

Theorem 5.6 Suppose that 7 G Mf , p > ~ and that ijj G W%, K, % G 
and jB Tj ^ > 0. Then the Lichnerowicz equation 115. 6\) admits a positive 
solution if > 0, <p G W%, in the following cases. 

(1) (7, VO zn i/ie positive Yamabe-scalar field class and Ayx^ ^ ; or 

(7, VO ^ e zero Yamabe-scalar field class and inf B T ^ > 0. 

Proof. If (7, %p) is in the positive or zero Yamabe-scalar field class a 
constant supersolution can be constructed directly as follows. First note 
that, as indicated in Theorem 15.5} we may assume that IZ^^ > 0. This 
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may require that we make a preliminary conformal transformation of our 
initial data, and find a new solution of the conformally formulated momentum 
constraint equation (see [8]). We consider the function of the single variable 
y defined by 

F(y) = B Tii! y^ + TZ^y^ - A^ n , (5.10) 
where / denotes the mean value of a function / on (E, 7) : 

vrikr) L fivoL <- 

3n-2 

Note that by setting y{x) = (fi(x) we see that f(x,4>) = y 4 <™- 2 ) F(t/), pro- 
vided that we do not replace the coefficients by their average values. 

Under the stated hypothesis one easily sees that F{y) is monotonically 
increasing on M + and has exactly one positive root. We let y$ = tp^ denote 
this root, so that F(y ) = 0. Now consider the linear equation 

3n-2 n + 2 

A y v = K^(p - Ayjt^o n - 2 + B T ^- 2 . (5.11) 

By our choice of (po the right hand side of this equation has mean value 
zero and is therefore orthogonal to the constants. Thus we may consider 
the function v G W%, with mean value zero on S, which solves (15. lip . The 
function 

<p + = <p + v — inf v > ip , A 7 </? + = A 7 v, 
is a supersolution if 1Z 1 ^ > because it holds that: 

3n-2 3n — 2 n + 2 n+2 

A^ + -f(-,<p + ) = K^(<po-<P+)-A ljk>7V ((p - 2 -<p + n ~ 2 )+B T ^- 2 -r+ 2 ). 

Hence if TZ 1 ^ > 0, then since A-g ,n — an d > 0, we have 

A 7 ^ + -/(-,^ + )<0 

because ip + > ipQ. In the case of the zero Yamabe-scalar field class the same 
type of argument holds, but we must use, in addition, the hypothesis B T ^ ^ 
to insure that ip > 0. 

In order to find a positive subsolution first note that any number £ < 1 
such that 

sup s (7?. 7iV , + B T ,f) 
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is a constant subsolution. It is positive only if inf s A 7 ^ n > 0. 

One may relax this hypothesis on A^ % n by instead constructing a non 
constant subsolution using the conformal invariance of the Lichnerowicz equa- 
tion [61 [TTJ [IH [8] . We ignore for the purposes of this theorem the connection 
between the coefficient A^ % n and the tensor K arising from the solution to 
the conformally formulated momentum constraint equation (15.51) . In order to 
state the conformal invariance properly one must, in addition to conformally 
rescaling the coefficients, solve (15. 5p with the conformally transformed data 
before posing the conformally transformed Lichnerowicz equation. We refer 
the reader to Proposition 2 of [8] for details. Since we only are concerned 
here with the Lichnerowicz equation we may state the conformal invariance 
at follows: 

3n— 2 n+2 

A 7 y? - TZ^^if + A^^ip "~ 2 - i3 TiV ,^™-2 = 

n + 2 / , i i .i , 3n-2 , , n+2 \ 

6—i (A 7 y - k' 7 ^' + A'^^y— - B' T ^'—z) , 

with 

y = ^ 7 , y/ = rV, A^ = A^ R J~^ and B' Tii) = B T ^. 
Now suppose A^ )7r > 0, A^^ ^ 0. We set 

k = TZj^ + XB T ^ 

with A = in the positive Yamabe-scalar field case, while we take A > 
(infs fv^) -1 in the case 1Z 1 ^ = 0. Then there exists a 9 > 0, 9 G W% , such 
that 

A,9 -k9 = -A^. 
Then A 7 # - K^9 = -9^11'^ implies 

K'^ = 0-& (A^ + (TZ^ - k)9) . 

The "primed Lichnerowicz equation" then admits the positive constant sub- 
solution £ if 

, , 3n — 1 , n+2 

or, equivalently, 

9~^{-(\k,* + ^)8}t + A 7A J-^r^ - B T ^ > 0. 
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Any number I such that 



n — 2 -1 n — Z 

£ < minfinf X~9 , inf 9^=^), 
s s 

is a positive subsolution of this transformed equation, and this shows that 
Q~ x t is a positive subsolution of the original equation. ■ 

An analogous method allows for the construction of sub and superso- 
lutions, in the negative Yamabe-scalar field class. We refer the interested 
reader to [3] and [16] for the vacuum case and [7j for the Einstein-scalar field 
system. 

Remark 5.7 The existence and non existence results presented here cover 
all cases when B T ^ > is a constant^ since then either inf^ B T ^ > or 
B T ^ = 0, and one is therefore studying the equation 

A</? > [or < 0] and ^ 0, 

which has no solution on a compact manifold. 



5 This is the case when the scalar field reduces to a cosmological constant. 
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